Spin waves in a Bose Ferromagnet 
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It is shown that the ferromagnetic transition takes place always above Bose-Einstein condensation 
in ferromagnetically coupled spinor Bose gases. We describe the Bose ferromagnet within Ginzburg- 
Landau theory by a "two-fluid" model below Bose-Einstein condensation. Both the Bose condensate 
and the normal phase are spontaneously magnetized. As a main result we show that spin waves 
in the two fluids are coupled together so as to produce only one mixed spin-wave mode in the 
coexisting state. The long wavelength spectrum is quadratic in the wave vector k, consistent with 
usual ferromagnetism theory, and the spin-wave stiffness coefficient c s includes contributions from 
both the two phases, implying the "two-fluid" feature of the system. c s can show a sharp bend at 
the Bose-Einstein condensation temperature. 
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Ferromagnetism belongs to the oldest phenomena in 
condensed matters and it is attracting continuous re- 
search interest till today .1], sparked by discoveries of new 
materials or new phenomena in this field. Most ferromag- 
nets being studied so far are composed of fermionic parti- 
cles. The recent experimental success of optical confining 
ultracold atomic Bose gases 0,0] provides possibilities to 
realize a new kind of ferromagnet, the Bose ferromagnet. 
A promising example might be the gas of F = 1 87 Rb 
atoms. It has been predicted in theory^ and confirmed 
by experiments^^ that the hyperfine spin-spin interac- 
tion in this Bose gas is ferromagnetic. 

In dilute atomic gases, interatomic forces are rather 
weak, with the effective s-wave scattering length being 
typically of the order lOOag where as is the Bohr radius. 
The spin-dependent interaction is even 1 or 2 orders of 
magnitude smaller 0. One may question whether the fer- 
romagnetic (FM) transition induced by such a weak FM 
coupling could be observed in experiments. Recently, 
it was already shown that the FM transition in Bose 
gases takes place always above Bose-Einstein condensa- 
tion (BEC), regardless of the value of the FM coupling^. 
It means that Bose gases can exhibit ferromagnetism at 
relatively high temperatures in comparison to the energy 
scale of the FM coupling. Below BEC, the FM Bose 
gas becomes a "two-fluid" system: the polarized Bose 
condensate coexisting with the magnetized normal gas. 
This is a unique feature in the Bose ferromagnet. 

Accompanying the FM transition, spin waves appear 
in the ferromagnet as a Goldstone mode. In conven- 
tional ferromagnets, both insulating and metallic, the 
spin-wave excitations are gapless at wave vector k = 
and the long wavelength dispersion relation is quadratic 
in k, u) s = c s k 2 , with the spin- wave stiffness coefficient 
c s proportional to the strength of the FM coupling^]. It 
is a rather interesting problem how the spin wave mani- 
fests itself in the Bose ferromagnet, especially considering 
the "two-fluid" feature of the system. In this letter, we 
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FIG. 1: Schematic Phase diagram of the Bose ferromagnet 
(left) and the SC ferromagnet ZrZn2 (right). In the left fig- 
ure, Tf is the FM transition temperature. T c and denote 
the BEC temperatures for the Bose gas with and without FM 
couplings (denoted as 7) respectively. In the right figure, Tf 
and T s are the FM and SC transition temperatures, respec- 
tively. The pressure P is supposed to modify the effective 
spin-exchanges between electrons in SC ferromagnets. 

shall examine the phase diagram and spin waves in the 
Bose ferromagnet phenomenologically and show how spin 
waves in the two fluids couple together. 

Spin waves in spinor Bose-Einstein condensates were 
intensively studied in theory]^, llCl 1 1 1| , motivated by 
experiments on spinor BEC|2(. In their pioneering pa- 
pers, Ho0 and Ohmi and Machida^ calculated low- 
lying collective modes both in fcrro- and in antiferro- 
magnetic condensates. For the FM case, a Bogliubov 
spectrum for the density fluctuation and a fc 2 -formed dis- 
persion for the spin fluctuation were derived, 



= V e l + 2 (-9o + 9s)n ek, hw Q = e k = -^—, (1) 

coinciding with available theories of BEC and ferromag- 
netism. Here go and g s denote the spin-independent and 
dependent interactions, respectively. A surprise in their 
results is that the spin-wave stiffness, c s = ft 2 /(2m), does 
not depend on the FM interaction. The same results 
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were obtained by Ueda through a many-body mean-field 
approach |10|. Also many attempts were made to gener- 
alize these results to finite temperatures 0. However, 
all these works are dedicated to only the Bose condensed 
phase, without taking into account the magnetized nor- 
mal phase. 

Analogous with the Bose ferromagnet, coexistence of 
ferromagnetism and superconductivity was discovered in 
various solid state materials 0, such as UGe2, ZrZn2 
and URhGe. It is observed that these materials first 
undergo a FM transition, then go into the supercon- 
ducting (SC) state. Their temperature-pressure phase 
diagram seems quite similar to the temperature-FM cou- 
pling phase diagram of the Bose ferromagnet 0, as shown 
in Figure 1. In such SC ferromagnets, both ferromag- 
netism and superconductivity are believed due to the 
same band electrons and the Cooper pairs are most likely 
in the triplet configuration^^ . Thus the system behaves 
like a spin-1 Bose gas to some extent. 

Phcnomcnological theory is proved to be a powerful 
tool to investigate behaviors of superconductors and su- 
perfluids near the transition temperature. It has already 
been successfully employed to describe various phase 
transitions |T^ and to explain the phase diagram [f! Il5| 
of the above-mentioned SC ferromagnets. and also has 
ever been applied to spinor Bose gases |lall7t . Although 
phenomenological theory only produces qualitative re- 
sults, it leads to deep insight into the nature of the asso- 
ciated phenomenon. 

We start with an appropriate Ginzburg-Landau (GL) 
free energy density functional allowing for the coexistence 
of BEC and ferromagnetism. It consists of three different 
parts: f t = f m + f b + f c , with 

f m = c|VM| 2 + ^|M| 2 + ^|M| 4 , (2a) 

f h = ^V*t. V * + a|*| 2 + §|*| 4 
2m 2 

+^*K'F tri ■ F^ 7 /* y * 7 , (2b) 

U = -gM ■ *;F ff7 * 7 , (2c) 

where repeated sub-indices represent summation taken 
over all the hyperfine states. This model is based on the 
two- fluid description of the Bose ferromagnet. / m (M) 
describes the ferromagnetic normal phase, which is the 
standard Landau free energy density for the second-order 
FM transition, having the magnetization density M as 
the order parameter 18] . /(,(*''', H*) describes the Bose 
condensed phase. \Pt = (v|> T )* is the complex order pa- 
rameter of the spinor condensate, m is the mass of the 
particle, c, a, 6, a, (3q and j3 § are phenomenological pa- 
rameters, c is the spin- wave stiffness in the FM normal 
gas. a and a depend on the temperature, a = a'(T — Tf) 
and a = a'(T - T c °) where T f and T c ° are the FM transi- 
tion and BEC temperatures in the decoupled case, respec- 



tively, a', a', b and f3o are positive constants and inde- 
pendent of the temperature. The (3 S term comes from the 
spin-dependent interaction and thus has SO(3) symme- 
try. F = (F x , F y ,F z ) are Pauli matrices for the hyperfine 
spin. (3 S < in a Bose ferromagnet. f c describes the cou- 
pling between the two phases, with the coupling constant 
.9 > 0. 

The full FM order parameter in the normal gas takes 
the form M = (M) + <5M. Its average value is cho- 
sen to be along the z direction for convenience, (M) = 
(0,0, M ). SM = (SM x ,SM y ,8M z ) represents spin fluc- 
tuations. The full order parameter for the condensate 
is written as * t = + <5**, $* + 5%, + 5V*_ 1 ). 
Q a = and = no is just the density of con- 

densed bosons. For a homogeneous system, both M and 
Q a are spatially independent, so VMo = = 0. 

First, we examine the phase diagram predicted from 
a microscopic model in Ref. 0- Neglecting all the fluc- 
tuations and minimizing the total free energy density ft 
with respect to , one gets the stable solution for the 
condensate. It is found that only spin-1 bosons condense, 

Po + p s V a' J 

and |$o| 2 = l^-il 2 = 0. These results agree with the 
prediction of Ho|g| and Ohmi and Machida[s|, and also 
the microscopic theory [?]. Obviously, the BEC critical 
temperature is enhanced by M , the magnetization in 
the normal phase, 

T c = T c ° + ^-M . (4) 
a' 

At T = T c , the order parameter of the condensate is zero, 
and we can derive the value of Mq by minimizing / m (Mo) 
with respect to Mo, 

M = \J y [Tf — T c ] . (5) 

Combining Eq. @ and JSJ, one has ST C — 

^/T*(ST f -ST c ), with T* = (g/a'fa'/b, 6T C =T C - T° 
and STf — Tf — T° . At very small ferromagnetic cou- 
plings, / — > 0, both 5T C and STf tend to zero with I. 
Provided that 5T f « T*[l||, we have 

iST c )*6T f (l-¥±y (6) 

Suppose STf depends linearly on / when I << l[l9|. 
STf = CI, the /-dependence of ST C is given by 

= CI (l-?Pj. (7) 

Besides a linear term, ST C also depends on / quadrati- 
cally. So far, the phase diagram of Bose ferromagnets is 
roughly reproduced. Similar treatment had been carried 
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out to explain the T — P phase diagram of SC ferromag- 
nets previously [Til Il5| . 

It was argued that BEC itself cause a FM transition 
in spinor bosons |io|. If one considers only the free en- 
ergy density of the Bose condensate, the /3= term alone 
indeed can break the 5*0(3) symmetry 0, 0|. However, 
since FM transition occurs already above BEC, the f5 s 
term no longer plays a critical role. It only enhances the 
condensate fraction, but does not result in another FM 
transition. This point is confirmed by the microscopic 
theory .7], which shows that both the magnetization M 
and the derivatives {dM/dT) seem smooth at T c . 

To proceed, we discuss spin waves in the Bose ferro- 
magnet. Within the GL theory, the spin wave is inter- 
preted as spin fluctuations around the average value of 
the order parameter. To linear order in <5M and 5^ a: the 
fluctuations are described by // = + /jf + with 



f 1 = 

H = 
fi = 



cVSM+VSM- + - (a + bM$) 6M+5M- , (8a) 
ft 2 

— v<y** v**„ + [a + (A, + a)*! 2 ]^;^ 

-2/3 s $ 1 2 ^*_ 1 ,5*_ 1 , 



-—g$i(6%5M+ + 5V 5M- 



(8b) 
(8c) 
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Here SM + = SM X + iSMy and SM + = SM X 
scribe the transverse spin fluctuations (i.e. the magnons) 
in the normal gas. The longitudinal spin fluctuation is 
neglected because it is a gapped mode in the FM phase 
and can hardly be excited at low energy Q. As pointed 
out in Refs. H S^>i, S^/ , and 5^_i represent the 
density, spin and "quadrupolar" spin fluctuations in the 
condensate, respectively, since Sn = v / no((5 , I'i + c^*), 
5M- = 0^*o and SAP = 2y/r^8^*_ x . Equations © 
can be simplified further making use of the conditions 
under which the second-order phase transitions occur, 
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0, 



- 5 M + a+(A,+/3 s )$i 2 = 0. 



(9) 
(10) 



In order to derive dispersion relations of the col- 
lective modes, Equations (JSJ should be expressed in 
the momentum space by performing Fourier trans- 
forms, <5* CT (r) = ^k<^( k ) ex P(* k ' r )> SM-(r) = 
2-Zk (k)exp(ik ■ r). For the density fluctuation of the 
condensate, we obtain 



/'(a*i,<y* a ) = $^(a*;(k) ^r(-k)) 



(ii) 



e fe +/3$i 2 / U*I(-k) 



5*1 (k) 



where j3 = (5q + j3 s and means the summation is 
taken only over k > 0. The frequency of S^i is given by: 



hux = \d + 2/3$^e fc , 



(12) 



which is of the Bogliubov from[^, Notet here (3 is 
a phenomenological parameter, not the real interaction 
strength. 

For the "quadrupolar" spin fluctuation, we have 
/ J (¥* !,*_!) = V<5*l 1 (k)( efe +2.gAf -2/3 s $ 1 2 )^_ 1 (k). 
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(13) 

Owing to the magnetization in the normal gas, the gap in 
the frequency of 6^f-i becomes larger. So basically, this 
kind of excitation is negligible in the Bose ferromagnets. 

Our main purpose is to derive spin waves (magnons) 
in the coexisting state of the Bose ferromagnet. It can 
be seen from Eqs. that the spin waves in the thermal 
could and in the condensate are coupled together, 

f{5M + ,5M^5%, 5* ) = (*M+(k) 5*5(k)) x 




2M 




ek + gMoJ V 5 *o(k) 



(.14) 



This equation indicates that the spin fluctuations in each 
phase become gapped solo (see the diagonal matrix ele- 
ments), due to the coexistence of the two phases. It seems 
violating the Goldstone theorem. To check this, we cal- 
culate the spectrum of the coupled spin mode, loq, which 
is determined by 



ck 2 



2M 



e k + gM + huj 



0. 



(15) 



There should be two solutions to Equation (|15J) . but only 
one solution is positive, 



Hujq 



(Ae) 



4{ck 2 e k + ^e k +gM ck 2 



where Ae = ck 2 



2M„ 



gM 



(16) 

It implies that 



only one mixed spin-wave mode exists inside the Bose 
ferromagnet although the system is separated into "two 
fluids" , different from the case of sound waves. This is 
physically reasonable, since spin wave is associated with 
the Goldstone mode and the FM transition happens only 
once. This mode mixes spin waves in the "two fluids" . 

Next let us look at the long wavelength dispersion of 
the spin wave. At very small k, the spectrum reduces to 



Hujq 



gM ck 2 



g±L_ 

2M 



Eft 



gM 



g<s>i 2 

2M 



(17) 
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Once again, we obtain the fc 2 -formed dispersion relation, 
and the Goldstone mode is recovered. We note that the 
coupling constant g is cancelled in the long wavelength 
spectrum. Since this result holds only in the vicinity of 
the BEC temperature, where the BEC order parameter is 
very small, $f << Mo, so the spectrum can be simplified 
further to 

1 $1 & 

fruin 1=3 C s k , with C s — C + —-t-tt;- — . (18) 

' 2 Ml 2m K ' 

c s is the spin-wave stiffness of the mixed spin- wave mode 
and it contains contributions from both the two fluids. 
But c s is not the simple summation of the spin-wave stiff- 
ness in the two phases. A weight factor defined by | jfe 
is applied to the condensate spin-wave stiffness. Usually 
c changes smoothly with temperature under Tf. Owing 
to the spin-wave stiffness of the condensate, there should 
be a sharp bend in c s at T c , which could be measured 
experimentally. 

In 2002, the JILA group observed spin-state segrega- 
tion in an ultracold, noncondensed Rb gas with two sub- 
spin states [2lJ. This effect can be interpreted as spin 
wave in a pseudo-spin-1/2 Bose gas 22]. In a more re- 
cent experiment, similar phenomenon was studied in a 
partially Bose-condensed gas^. Coherent coupling be- 
tween the normal and superfluid components was ob- 
served. Our work is relevant to the latter casejl^. As we 
predict above, the spin waves in two components should 
be coupled into one mixed mode (the coherent mode). 

For triplet superconductors, the order parameter of the 
SC phase is in a 2 x 2 matrix form, A(k) = i[d(k) • a]a y . 
The product A(k)A(k)t = |d(k)| 2 <7 -H[d(k) x d*(k)]-a, 
where a = (a x , a y , u z ) is the Pauli matrix and m(k) oc 
i[d(k) x d*(k)] acts as magnetic moment of Cooper 
pairs. It is suggested that the strong internal field in 
SC ferromagnets stabilize a nonunitary triplet state|l3|. 
which means that Cooper pairs carry a net magnetization 
m/0. This case is equivalent to the polarization of the 
FM spinor Bose condensate, and thus the system consists 
of two coexisting ferromagnetic phases. Seeing the simi- 
larity between the Bose and SC ferromagnets, our results 
for spinor bosons help to understand spin waves in SC 
ferromagnets. This problem has drawn current research 
interest |25J. The present phenomenological approach can 
be extended to SC ferromagnets. 

In summary, we have shown that the Bose ferromag- 
net consists of "two fluids" at low temperatures, namely, 
the magnetized normal and Bose-condensed phases. Nev- 
ertheless, only one mixed spin-wave mode exists inside. 
The spin-wave spectrum has the same momentum de- 
pendence as that in a conventional ferromagnet at long 
wavelengths, u> — c s k 2 . However, the spin- wave stiffness 
c s contains contributions from the two different phases. 
c s is a function of the condensate density no below BEC, 
so it can exhibit a sharp bend at the BEC temperature, 
which could be examined experimentally. 
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